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SUMMARY 



New error bounds are developed for the Liouville- 
Green approximation to the solution of an important class 
of differential equations arising in military operations 
research {specifically, variable-coefficient Lanchester- 
type equations of modern warfare for combat between two 
homogeneous forces) . In contrast to many previous re- 
sults, our error bounds apply to initial-value problems 
and are expressed in terms of initial conditions. Previous 
error bounds for boundary-value problems are sharpened as 
a consequence of our development of these new error bounds 
for initial-value problems. Finally, applications are 
made to some important specific models of combat between 
two homogeneous forces with time-dependent attrition- 
rate coefficients. 



1. Introduction 



OLVER [12] has developed error bounds for the so-called 
LIOUVILLE-GREEN^^ (LG) approximation [5, 10] to the solution of 
the differential equation 

= J(t)x. (1) 

His results extended earlier work by BLUMENTHAL [1] . The LG 
approximation is of particular importance in applied mathematics 
because it involves only elementary functions. As OLVER [12] 
has pointed out, however, the development of strict upper bounds 
for the errors in the approximate solutions has received rela- 
tively little attention. Moreover, the error bounds that have 
been previously developed [12, 13] are for boundary- value prob- 
lems, since most problems of interest in mathematical physics 
are boundary-value problems . These error bounds do not apply 
to initial-value problems. 

Thus, the purpose of this paper is to develop error bounds 
for the LG approximation to initial-value problems. Furthermore, 
in recasting OLVER' s results in a form suitable for initial-value 
problems , we have been able to sharpen his bounds for boundary- 
value problems. Although we will develop our results within 
the context of a specific problem in operations research, they 
are clearly applicable to the general second order initial-value 
problem, and it is a straightforward task to so recast them. 

^^Also called the WKB approximation [13]. 



Probably the most widely used (at least in the United 

States) deterministic differential-equation model in operations 

2 ) 

research [4,15-17] are LANCHESTER-type equations of warfare , 
which are named for the pioneering 1914 work of F. W. LANCHESTER 
[9]. In this paper we consider such a linear, variable-coeffi- 
cient differential-equation model for combat between two homogene 
ous forces. This model yields an X force-level equation equiva 
lent to (1) . Unfortunately, for even the simplest time-varying 
attrition-rate coefficients of interest, the solutions cannot 
in general be expressed in terms of either elementary functions 
or tabulated higher transcendental functions [17] . It is there- 
fore natural to seek an approximation in terms of elementary 
functions. Moreover, one is always interested in a simple 
a priori estimate for the error in the approximate solution. 

In this paper we develop an error bound that is both realistic 
and easy to evaluate. This bound is then computed for some 
particular attrition-rate coefficients of interest. 

2 . Variable-Coefficient LANCHESTER-Type Equations of Modern 
Warfare 

We consider the following variable-coefficient LANCHESTER- 
type equations of modern warfare^^ for combat between two 
homogeneous forces 



2 ) 

We will refer to any differential-equation model of combat as 
being LANCHESTER-type equations. The state variables are 
typically the force levels of the various weapon system types. 

3) 

The term "of modern warfare" denotes that we are considering 
linear differential equations. There are other nonlinear types 
of LANCHESTER equations [4,16]. 
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dx/dt = -a{t)y, dy/dt = -b(t)x, (2) 

with initial conditions 

x(t=0) = Xq, y(t=0) = y^ 

where t = 0 denotes the time at which the battle begins, x(t) 
and y(t) denote the numbers of X and Y at time t, and 
a(t) and b(t) denote time-dependent LANCHESTER attrition- 
rate coefficients. These equations (2) have been hypothesized 
to model combat in which both sides use aimed fire and target 
acquisition times are independent of the numbers of firers and 
targets [2,19]. The attrition-rate coefficients represent the 
effectiveness of each side's fire (i.e. its firepower). Temporal 
variations in a side's fire effectiveness (caused by, for example, 
changes in force separation, combatant postures, target acquisi- 
tion rates, etc.) are modelled by the time-dependent attrition- 
rate coefficients. Further discussions of the physical assump- 
tions hypothesized to yield (2) , estimation of the attrition- 
rate coefficients, and the importance of (2) to military opera- 
tions research are found in [16,17]. 

We assume that a(t) and b(t) are positive and twice 
differentiable for t^ < t < +<» with t^ < 0. We also assume 
that a(t), b (t) eL(tQ,T) for any finite T. We further take 
a(t) and b(t) to be given in the form a(t) = k g(t), 
b(t) = k h(t), where k , k are positive constants chosen so 

O ^ D 

that a(t)/b(t) = k /k. when g(t) = h(t) for all t. We 

3i O 

introduce (see [17]) the intensity of combat I(t) and the 
relative fire effectiveness R(t) defined by 

I(t) = /a(t)b(t) , and R(t) = a(t)/b(t). (3) 
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and 



We similarly introduce the combat-intensity parameter 
the relative-f ire-effectiveness parameter X^^ defined by 




and 




(4) 



A large class of tactical situations of interest can be 



modelled with the following general power attrition-rate 



coefficients [17] 



a(t) = k (t+C)^ 



and b(t) = k, (t+C+A) 

D 



V 



(5) 



where A,C >0. We will call A the offset parameter , since it 
allows us to model (with p,v>0) battles between weapon systems 
with different maximum effective ranges. We will call C the 
starting parameter ^ since it allows us to model (again, with 
li,v>0) battles that begin within the minimum of the maximum 
effective ranges of the two systems. The offset and starting 
parameters are related to various physical quantities in [17] . 

We observe that t^ = -C. Also, a(t) , b(t) £L(tQ,T) implies 
that y , V > -1 . 



right-most finite singularity of the X force-level equation. 
Furthermore, we set = 0 if there are no finite singularities. 



From (2) we may obtain the X force-level equation 




2 



( 6 ) 



with initial conditions 



x(t=0) = X 



0 



and {[l/a(t) ]dx/dt)^^Q = -y^. 



X Y X 

We may consider that t^ = max(tQ,tQ), where t^ denotes the 
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3 . LIQUVILLE-GREEN Approximation to LANCHESTER-Type Equations 
of Modern Warfare 

Let 



and denote T(t=0) as t^. Then Tq - 0 

a(t), b(t) eL(tQ,T) it follows that t = T(t) is well defined 
by the CAUCHY-SCHWARZ inequality for integrals. The transforma- 
tion is also easily seen to be invertable. We observe that 
T - Tq is related to the average intensity of combat I (t) by 



We may call t - "the elapsed normalized battle time," since 
the transformation (7) reparameterizes the battle's time scale 
in terms of elapsed time of battle and average combat intensity 
as (8) shows us. The siabstitution (7) transforms (6) into 




(7) 




( 8 ) 




(9) 



with initial conditions 



x(t=Tq) = Xq, and { [1/R^^^ (t) ] dx/dx 



0 



Remark 1: It is easily shown (see [17]) that (9) may be 



transformed into a linear second order differential equation 



with constant coefficients if and only if 



1 ^ 
I(t) dt 



^2.n R{t) 



constant. 
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Let Rq denote R(t=0) . The si±)s titution 

x(t) = X(t) [R(t) (10) 

transforms (9) into LIOUVILLE's normal form (see INCE [7] or 
KAMKE ( 8 ] ) 

1^- {1+F(t) }X = 0, (11) 

with initial conditions 

X(t=Tq) = Xq, and dX/dT(T=XQ) = -Yq/R^-XqEq, 

where 

F(t) =P"(t)/P(t), P(t) = [R(t)]“^/^, (12) 



41 (t) dt 

£q denotes e(t=0), and P' (t) denotes dP/dx. 

Writing (11) as d^X/dx^ - X = F(x)X, we may use variation 
of parameters to obtain the solution to (11) as 



X(x) = Xq cosh (t-Xq) - (Yq/^+x^Eq) sinh(x-XQ) 



[ F (a) si 
x« 



inh (x-o)X(a)da. (14) 



.4) 



If one drops the integral term in (14) , one obtains the 
LIOUVILLE-GREEN approximation X(x) 



4) . . 

Heuristically , if the appropriate fractional power of the relati 
effectiveness R(t) is "slowly varying," then by (12) we would 
expect that |f(x) |<<1 so that the integral term in (14) is 
"negligible." Theorem 1 gives us bounds on how "negligible" 
this term is. 
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( 15 ) 



X(t) = Xq cosh (t-Tq) - (y^y^+x^eQ) sinh (t-Tq) , 
which in terms of the original independent variable x reads 



X 



(t) = {Xq cosh (t-Tq) - (yQ/R^+x^eQ) sinh(T-TQ) } . (16) 



We observe from (14) that F(t) >0 for all ^ implies that 

as long as x(t) >0 we have x(t) > x(t) . A similar statement 
holds for F(t) f 0. As we shall see below, such cases in which 
F(t) is always >0 or 5 0 are readily encountered in applica- 
tions . 

4 . Error Bounds for the LIOUVILLB-GREEN Approximation 
The main result of this paper is Theorem 1. 

Theorem 1: Error bounds for the LIOUVILLE-GREEN approximation 

are given by 



|x (t) -x(t) I 5 Xq Kj e (t) < Xg Ky e (t) , 



(17) 



where 



= 2{(l+|eg|) +5^'^) , 



(18) 



V, 



J = I for 1-3 ^v^<£q and = 1 + ^ (19) 



y y 

J = II for -1 2 . < e„ < 1 2 and * 2, (20) 



Xq 0 0 



X(, 0 



II 



J = III for and ^ ^ = 1 - e„ - |2 0 , (21) 



and 
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e(t) 



( 22 ) 



{exp(^[ I F (a) I da) -llsinh (t-t ) . 

u Rq -I 2 J Tq 0 

The sign of the error is determined by the sign of F(t) . As 
long as x(t) >0, it follows that 

F(t) > 0 for all x ^ 0 implies that x(t) ^ x(t) , 

with the last inequality being reversed when F(t) 50 . 

Proof: Theorem 1 readily follows from Lemma 2, which is proven below. 

5 . Development of Error Bounds 

Consider the following fundamental system of solutions 
{X^,X 2 > to (11) 

X,^(t) = (l+hj^(T) }exp( (-1 )^^”^^t -Tq) ] for k = l,2, (23) 

where h, (x) for k = 1,2 is to be chosen so that 
k 

hk(x=To) = 0, and dhj^/dx (t=Tq) =0. (24) 

It follows that the solution to (11) may be expressed as 

X(t) = ^{(xQ(l-eQ)-yQV^)e ^ (l+h^(T)] 

- (T-T ) 

+ (Xq (I+Cq) + yQ/^) e [l+h 2 (T))}, (25) 

so that 

X(t)-X(t) = ±{[xQ(l-£Q)-yQ/R^]e ^hj^(x) 

-(t-t ) 

+ [Xq (l+£Q)+yQ»^]e ^ h2(T)}. (26) 
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bstituting (23) into (11), we find that ^ 

St satisfy 



d^h 



k k '^’’k 

3t^- 



( 27 ) 



th initial conditions (24) . 

We may consider hj^(x) for k = 1,2 to be an error term 
r the LIOUVILLE-GREEN approximation. We next develop a bound 
its magnitude, which sharpens earlier results by OLVER [12]. 



mma 1 ; A bound on the magnitude of ft>r k = 1,2 is 

yen bv 

t) |exp ((-1)^^"^^t-Tq) ] < 2(exp(^ [ 1 F (a) | da) -1} sinh ( t-t^) . (28) 

J 'T 
^0 



oof; For notational convenience, we develop the bounds for 
(t) separately. Transposing the right-most term on the left-hand side of 
7) for k = 1, treating { 1+h^ (t) }f( x) as a "forcing term," and inte- 
iting twice; we obtain the following VOLTERRA integral equation after a 
rther integration by parts 






h^(T) =- 



}f(o) {l+h^(o) }do. 



(29) 



Lving (29) in the usual manner by successive approximations, we obtain 

00 

h, (t) = y T (t) , (30) 

1 n 

n=l 

ire Tq ( T ) = 1 and for n > 1 



C (T) =^f }F(o)T ,(o)do. 

n 2 I n-1 



(31) 
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Observing that { 1-e 



we find that 



^ ^ }51-e for 0<T^Ha<T/ 



|t (T)| <ii;2_ L I 



F(a) |da , 



with equality holding for ^ straightforward inductive argument 



along the usual lines now yields 



n 



1t^(t)| ^ 



2 n! 



lF(a)ldo> . 



(32) 



0 



One step in the inductive proof of (32) is deserving of further elaboration, 
however. From (31) and (32) we obtain 

-2(a-T^)s ^ rO 



|t (T) 1 <— — 
n+1 , 



{1-e 



2(a- 



, -2ta-T 

nF(a)K |f(u) ]du> da. (33 



The inductive proof of (32) is completed by combining (33) with the observe- 

-2 (a-T^)l 



X.- 4-u ^ 2 (a-T) . 

tion that il-e i 



-2(t-To> 

e for 0<T^5a5x. Our 



sharpening of OLVER's results is due to this observation. The remaining 
steps in the proof of (28) for k = 1 follow along well-known lines [12,13] 
and will be omitted here. Similar arguments are used to prove (28) for 
k = 2. 



Remark 2: In our notation, OLVER's [12] corresponding 

error bound for k = 1 would read 



|h (t) I _< exp(^ / |f(o) |da)-l . 

i ^ Tq 

His corresponding error bound for k = 2 is not directly 



(34) 



comparable to our result here, since he does not take 
both errors zero at the same point. 



5 ) 



In our notation for the finite or infinite interval 
OLVER [12,13] takes h^(T=TQ) = dh^/dx(T=TQ) = 0 and 
h 2 (T=Tj^) = dh^/dT ( T=T ) =0 in contrast to (24). This is 
what makes his results unsuitable for initial-value problems. 



Remark 3 : Similar argiiments may be used to develop a 

bound on | h ’ ( t ) | . For our applications this result is 
not important. 



Lemma 2 : In terms of the transformed dependent variable X(t) / 

error bounds for the LIOUVILLE-GREEN approximation are given by 



|x(t)-X(t) I < Xq Kj E(t) < Xq Ky E (t) , 



Proof: From (26) and Lemma 1 we obtain 



(35) 



where and Kj are given by (18) through (21 ) , and 

E(t) = {exp(i| |F(a) |da)-l}sinh (t-Tq) . (36) 



|x(T)-X(T) I < t UQ(l-eQ)-y^/R^| + UQ(l+e^)+yQ»^| }e(t) 



(37) 



It follows that a rather loose error bound is given by 



|x(T)-X(T) I < XqK^E(t) . 
y y 

We observe that for -1 < 1 have 



Xq 0 -0 









(38) 



(39) 



so that (37) becomes | X(t) -X(t) | 5 2XqE(t) . Thus, (35) is proven for J = 

Yq 

For 1- — e , the error bound (37) becomes 

Xq u 0 



0^ jf 

|x(t)-X(t) I < 2X- (e- +— v^)E(t) , 

' ' 0 0 Xq 0 



(40) 



since 






( 41 ) 
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The error bound (40) may be sharpened, however, as follows. If F(o) >0 fo; 



0 < X 5 o 5 T, then as long as X(x) >0 we have X(x) > X(x) and g, (x) > 0, 
0 K 



where for k = 1,2 



= yhj^(x) exp[ (-I)^(x-Xq) ] . 



(42) 



It follows from (26) that for 1 — — /r^ — e. 

Xq 0 0 



0^ x(x) - x(x) < {xq(i+£q) ' 



(43) 



where 



E^(x) = {exp(j f F(o)da) - l}sinh(x-x^) . 

X. 



(44) 



Similarly, if F(o)^0 for x^5 0$x, then as long as X(x)>0 we have 



X(x)SX(x) and g (x)20, so that 



0> X(x) - X(x) > {xQ(l+eQ)+yQ.^}E_(x) , 



(45) 



where 



E (x) = {1-exp ( 



^ X, 



F(a)da) Isinh(x-XQ) . 



(46) 



It follows from (43) and (45) that 



|x(x)-X(x)| < x-{l+ e. +— »^}E(x) 

' ‘ 0 0 Xq 0 



(47) 



y y y 

Furthermore , since > 1 — implies that 2 (e^ + — »'^^) > 1 + e,, + — 

0 Xq 0 0 Xq 0 0 Xq 0 



the bound given by (47) is sharper than that given by (40). Thus, (35) is 
proven for J = I. The proof of (35) for J = III is similar to that for 



J = II. It is easily seen that K < K for x„,y > 0. 

J U 0 ■'o 
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6 . Examples 



We now compute theoretical error bounds for two special 
cases of the LIOUVILLE-GREEN approximation to the solution of 
(6) with the general power attrition-rate coefficients (5); 

(I) power attrition-rate coefficients with no offset (i.e. A=0), 

and (II) linear attrition-rate coefficients with positive offset 
(i.e. A > 0) . 

6.1. Power Attrition-Rate Coefficients with No Offset 
In this case we have 

a(t) = k^(t+C)^, and b(t) = kj^(t+C)^, (48) 

with C>0 and y,v>-l. The LIOUVILLE-GREEN approximation to 
the X force level is 

x(t) = (1+t/C) ^^^^{x^cosh (t-Tq) 

+ c''*)sinh(T-Tg)), {49) 

where 

T(t) = (l/6)Xj(t+C)'^, (50) 

and 

6 = (y+v+2)/2. (51) 

In preparation for estimating the error in the LIOUVILLE- 
GREEN approximation (49) by Theorem 1, we compute 

F(T) = (P“V) (3uyt4) ^ (32, 

166^t 
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Thus, F(t) > 0 for all t > Tq > 0 if and only if y > v. Hence, 
as noted at the end of Section 3, one frequently encounters in 
applications cases in which F(t) always has the same sign. 

For the error estimate (17) , we then have 

i f |F(o)|da = (53) 

Remark 4: The exact solution x(t) to (6) with attrition- 

rate coefficients (48) is given in TAYLOR and BROWN [17] (see 

also [16]). It may be expressed in terms of modified Bessel 

functions of the first kind of (for y,v>-l) fractional order, 

i.e. I for 0 < a < 1. Since few of the latter are tabulated 
a 

(i.e. tabulations only exist for a = ±1/4, ±1/3, ±1/2, ±2/3, 

± 3/4, and these do correspond to cases of interest) , TAYLOR 

and BROWN [17] suggested the use of new transcendents which 

•• 6 ) 

they called LANCHESTER-CLIFFORD-SCHLAFLI functions . 

6.2. Linear Attrition-Rate Coefficients with Positive 
Offset 

a(t) = k^(t+C), and b(t) = k, (t+C+A), (54) 

d JD 

with A,C> 0. The LIOUVILLE-GREEN approximation to the X forc« 
level is 



6 ) 



After earlier work by W. K. CLIFFORD [3] and L. SCHLAFLI [14] 
(see also [6,18] ) . 
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x(t) 



= [1 



(1+A/C) 
(l+A/(t+C) ) 



-] 



1/4 



{x^ cosh (t-t.) 



y - /XT 

■^0 R 



XqA/C 



yi+A/C 4XjC (1+A/C) 



3/2 



sinh(T-TQ)}, (55) 



where 



T(t) = Hn(ij;+/5p^) }, 



(56) 



and 



ip (t) = 1 + 2 (t+C)/A. 



(57) 



For estimating the error in the LIOUVILLE-GREEN approxi- 
mation (55), it is more convenient to express F(t) as defined 

by (12) in terms of the original independent variable t. We 
7) 

compute that 



A{12 (t+C) + 7A) 
16Xj(t+C) ^(t+C+A) ^ 



(58) 



It follows that F(t) >0 for all x > Tq > 0. Thus, 
previous example, F(x) always has the same sign, 
error estimate (17) , we then have 

0<j| F(a)da< min (m^ (t) ,m^ (t) ) , 

^0 



where 



as in the 
For the 



(59) 



m^ (t) 
m2 (t) 



3A 

16C^ (1+A/C) 

^ {4(l-q^) + 

32C^Xj 



{4(l-ql/2j 



■7A/T ,,4 , 

4C ^ ' 



^ 7A,^ 3/2, 1 

+ 4c(i-q 



(60) 

( 61 ) 



7 ) 



In general 
F(t) 



4b^(t) 



dr d 0 
dt dt 



b(t) - 



1 d „ 

T jf 
4 dt 
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and 



q(t) = 1/U + t/C) . 



(62) 



It may be shown that m^(t^) = m 2 (tj^) implies that m^(t) > m 2 (t) 
for all t >t^. The error control term estimate m^ (t) was 



Remark 5 : This case is of particular interest in military 

operations research, since it may be used to study combat 
between two weapon systems with different maximum effective 
ranges [16]. The exact solution x(t) to (6) with attrition- 
rate coefficients (54) is given in [17] (see also [16]). It 
apparently cannot be expressed in terms of previously "known" 
transcendents, since the X force- level equation in this case 
could not be found to correspond to any second order linear 
equation considered in [8] or [11]. 



7 . Final Remarks 

Although given within the context of a specific problem 
in operations research, the reader will have no trouble trans- 
lating the above results into those for the general second order 



developed for "small" t^^ , while m 2 (t) for "large" t^^ 



8) 



Letting 




we find that 
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Here we have used the bound J < 




9) 



In this case (to be contrasted with the previous one) , we have 



used the bound J s 



I 



^ {12 (s+C)+7A}d s 
0 (s+C)^ 
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initial-value problem. In the examples of Section 6 we saw 
that for two models of considerable importance in military opera- 
tions research both the LIOUVILLE-GREEN approximation and bounds 
on its error were simply expressed in terms of elementary func- 
tions. No previous application of the LIOUVILLE-GREEN approxi- 
mation has appeared in the operations research literature. In 
a subsequent paper we plan to present a numerical investigation 
of the accuracy (both numerical eveluation of the theoretical 
error bounds and a comparison with the exact solution when 
available) of the LIOUVILLE-GREEN approximation to the solution 
of LANCHESTER- type equations of modern warfare for combat between 
two homogeneous forces . 
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